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Abstract 



The present study deals with spatially homogeneous and anisotropic locally rotationally symmetric (LRS) 
Bianchi type I cosmological model with dominance of dark energy. To get the deterministic model of Uni- 
verse, we assume that the shear scalar (er) in the model is proportional to expansion scalar (0). This condition 
leads to A — B n , where A, B are metric potential and n is positive constant. It has been found that the 
anisotropic distribution of dark energy leads to the present accelerated expansion of Universe. The physical 
D ' behavior of the Universe has been discussed in detail. 
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, &* • 1 Introduction 

£NJ . The discovery of the accelerated mode of expansion of the Universe stands as a major breakthrough of the obser- 
vational cosmology. Survey of cosmological distant type la supernovae (SNe la; Riess et al 1998; Perlmutter et 
al 1999) indicated the presence of a new unaccounted-for Dark energy (DE) that opposes the self-attractions of 
matter and causes the expansion of Universe to accelerate. This acceleration is realized with negative pressure 

|p^ and positive energy density that violate the strong energy condition. This violation gives a reverse gravitational 

effect. Due to this effect, the Universe gets a jerk and the transition from the earlier deceleration phase to the 
recent acceleration phase takes place (Caldwell et al. 2006). The cause of this sudden transition and the source 
of accelerated expansion are still unknown. The state of the art in cosmology has led to the following present 
distribution of the energy densities of the Universe: 4% for baryonic matter, 23% for non baryonic dark matter 
and 73% so-called DE (Spergel et al. 2007). 



The isotropy of the cosmic microwave background (CMB) radiation, first seen by the cosmic background 
explorer (COBE) satellite (Smoot et al. 1992) and then reinforced by the Wilkinson Microwave Anisotropy 
Probe (WMAP) data (Hinshaw et al. 2003), together with the assumption that we are not in spacial position in 
the Universe, underlines the cosmological principles, according to which we live in a homogeneous and isotropic 
Universe described by a FRW line-element. Tiny deviation from perfect isotropy at the level of 10 -5 , have also 
been reported by Bennett et al (1996) and thereafter confirmed by high resolution WMAP data. The observed 
CMB anisotropy spectrum is in impressive agreement with the predictions of ACDM model. Koivisto and Mota 
(2008a, 2008b) proposed the mechanism of DE with anisotropic equation of state (EoS) parameter which is very 
attractive because cosmic anisotropy originates from the actual dominant component of the Universe and then 
could be directly tested, for example, by either observations of the magnitude and redshift of type la supernovae 
or cosmic parallax effects of the distance source. DE has been conventionally characterized by the equation of 
state (EoS) parameter u/ de ) = p( de )/p( de ) which is not necessarily constant. The simplest DE candidate is the 
vacuum energy (w = —1), which is argued to be equivalent to the cosmological constant (A) (Martins, 2002). 
However, it is well known, there are two difficulties arising from the cosmological constant scenario, namely the 
two famous cosmological constant problems - the fine tuning and the cosmic coincidence one. An alternative 
proposal is the concept of dynamical DE. Such a scenario is often realized by some scalar field mechanism and 
suggests that the energy form with negative pressure is provided by a scalar field evolving under a properly 



constructed potential. So far, a large class of scalar-field DE models have been studied, including quintessence 
viz a/ de ) > —1 (Steinhardt et al 1999), phantom viz u/ de ) < —1 (Caldwell 2002) and quintom (that can across 
from phantom region to quintessence region) . The quintom scenario of DE is designed to understand the nature 
of DE with w( rfe ) across -1 (Setare 2006). Recently, Cai et al (2010), Setare and Saridakis (2007, 2009b) have 
studied the DE models with EoS parameter across -1 which give a concrete theoretical justification for quintom 
paradigm. In addition, the other proposals on DE include interacting DE model (Setare 2007) and braneworld 
model (Setare and Saridakis 2009a) etc. By combining data from seven CMB experiments with large scale struc- 
ture data, the Hubble parameter measurement from the Hubble space-telescope and luminosity measurements of 
SN la, Melchiorri et al (2003) demonstrated the bound on w (de) to be -1.38 < uj {de) < -0.82 at 95% confidence 
level. 

The simplest of anisotropic models are Bianchi type-I homogeneous models whose spatial sections are fiat 
but the expansion or contraction rate are direction dependent. For studying the possible effects of anisotropy in 
the early Universe on present day observations many researchers (Huang 1990; Chimento et al. 1997; Lima and 
Troden 1996; Lima and Maia 1994; Pradhan and Singh 2004; Pradhan and Pandey 2006; Saha 2006a, 2006b) 
have investigated Bianchi type-I models from different point of view. The binary mixture of perfect fluid and 
DE has been studied for Bianchi type I (Saha 2005). Some Authors (Akarsu and Kilinc 2010; Yadav and Yadav 
2011; Yadav et al 2011b; Kumar and Yadav 2011; Amirhashchi et al 2011 and recently Yadav 2011) have studied 
anisotropic DE models with constant deceleration parameter (DP). 

In this paper, we considered minimally interacting perfect fluid and DE components with proportionality 
relation between shear scalar and expansion within the framework of LRS Bianchi-I space-time in general rel- 
ativity. The paper is organized as follows: In section 2, the models and field equations have been presented. 
Section 3 deals with the exact solutions of field equations and physical behavior of the model. The statefinder 
and distance modulus curves are described in section 4. Finally the results are discussed in section 5. 



2 The Metric and Field Equations 

We consider the LRS Bianchi type I metric of the form 

ds 2 = -dt 2 + A 2 dx 2 + B 2 (dy 2 + dz 2 ) , (1) 

where, A and B are functions of t only. This ensures that the model is spatially homogeneous. 

The Einstein's field equations in case of a mixture of perfect fluid and DE components, in the units 87rG = 
c = 1, read as 

R> J -±g i J R=-T< j , (2) 

where T*- = 1 ■ l + T ■ L is the overall energy momentum tensor with T^- and T i ■ as the energy momentum 
tensors of ordinary matter and DE, respectively. These are given by 

T ( ™ )1 = diag[- P {m \ P {m) ,p (m) ,p (m) ] , =<Ba0[-l,w< m \w( m >,w< wl V wO ( 3 ) 

and 

T { f l = diag[-p {de \p {de \p {de \p {de) ] = diag[-l,^ de \J de \J de ^p^ (4) 

where p^ m ' and p( m > are, respectively the energy density and pressure of the perfect fluid component or ordinary 
baryonic matter while u^ m ' = p\ m ) J p( m ) is its EoS parameter. Similarly, p( de ) and p( de ) are, respectively the 
energy density and pressure of the DE component while oj( de > = p( de ) /p( de ) is the corresponding EoS parameter. 
The Einstein's field equations @ for the line-element ([TJ lead to the following system of equations 

2- + i!J = -« (m v m) - u {de) P {de) , (5) 

B B l 

4 + 1 + 4^ = - w (™y™) - u;^p^ , (6) 

A B AB v ' 



|f + 2^§=P (m) +P (de) . (7) 

The Bianchi identity G ?j = yields 

p {m) + 3(1 + w (m) )p W i? + P (de) + 3(1 + uj ( - de '>)p ( - de ')H = 0, (8) 

with H being the mean Hubble parameter, which for LRS Bianchi I space-time can be defined as 

»-S-§(4+>§) 

Here, and in what follows, over-dots indicates differentiation with respect to t and a is the average scale factor 
of LRS Bianchi type I model: 

a={AB 2 )i. (10) 

The spatial volume (V) is given by 

V = a 3 = AB 2 . (11) 

The expansion scalar (0), shear scalar (a) and mean anisotropy parameter (A m ) are defined as 

e = 3H=^ + 2^, (12) 



^(xx-H. 



(13) 
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3 Solutions of the Field Equations 

In order to solve the field equations completely, firstly we assume that the perfect fluid and DE components 
interact minimally. Therefore, the energy momentum tensors of the two sources may be conserved separately. 
The energy conservation equation T *™ y = 0, of the perfect fluid leads to 

p (m) +3(l + u;( m V m )# = 0, (15) 

whereas the energy conservation equation T y = 0, of the DE component yields 

p {de) +3(1 W de) V (de) # = • (16) 

Following Akarsu and Kilinc (2010), we assume that the EoS parameter of the perfect fluid to be a constant, 
that is, 

- M = -^ = canst., (17) 

while uj( de *> has been allowed to be a function of time since the current cosmological data from SNIa, CMB and 
large scale structures mildly favor dynamically evolving DE crossing the phantom divide line (PDL) as discussed 
in Section 1. 

Finally, we constrain, the system of equation with proportionality relation between shear (a) and expansion 
(9). This condition leads to the following relation between the metric potentials 

A = B n (18) 

where n is positive constant. For anisotropic model n =/= 1. 
Equations ©, © and JTSJ lead to 

! + (n + l)|J = 0. (19) 



The solution of equation (|!20|) is given by 

B = (ht + ko)^, (20) 

Where fco and k\ are the constants of integration. 
From equations (TT8|) and (|20|) . we obtain 

A=(fcit + /*>)*. (21) 

The rate of expansion in the direction of x, y and z are given by 

i _ nfci 1 

^ i4"(»+2)(M + *0)' l ' 

H v = H *=B = (n + 2)(fc!t + A: ) ' (23) 

The mean Hubble's parameter (if), expansion scalar (0) and shear scalar (c) are given by 

H = o,, l l\,s , (24) 

3(fcii + fcoJ 

" fcl (25) 



(M + fco) 



9 (n-l) 2 fc 2 1 , s 

** = ,,„ ,, a,, ^ 2 - ( 2 6) 



3(n + 2) 2 (fcit + fco) 2 ' 
The spatial volume (V), mean anisotropy parameter (A m ) and DP (q) are found to be 



From equations (|23|) and (|2"o) . we obtain 



V = 


(M + feo) 


A m = 


2(n-l) 2 


(n + 2) 2 


"-( 


-1-1- 


d£ \ 


vW 


<7 


(n-l) 



(27) 
(28) 
(29) 



9 V3(n + 2)' (30) 

It is important to note here that the proportionality relation between shear and expansion leads to the positive 
deceleration parameter (q) with isotropic distribution of DE in LRS Bianchi -I space-time. Since we are looking 
for a model explaining an expanding Universe with acceleration, so, we assume the anisotropic distribution of 
DE to ensure the present acceleration of Universe. Thus equations (O, ([6]) and (fT6)) may be re- written as 

2^ + |- = -wWpW - {uj^ + S)p {de) , (31) 

B B z 

~A + I + A^ = - w{m)p{m) ~ {u{de) + l)p(de) ' (32) 



P 



f 3p (de) (l + ^ {de) )H + p {de) (5H x + 2 1 H y ) = . (33) 

The third term of equation (l33l) arises due to the deviation from a/ de ) while the first and second terms of 
equation (|33|) are deviation free part of T ■ \ According to equation (|33|) . the behavior of p^ de ^ is controlled 
by the deviation free part of of EoS parameter of DE but deviation will affect p( de ) indirectly, since as can be 
seen later, they affect the value of EoS parameter. But we are looking for physically viable models of Universe 
consistent with observations. Hence we constrained S(t) and j(t) by assuming the special dynamics which is 
consistent with (|33p . The dynamics of skewness parameter on x-axis (5) and y-axis or z-axis (7) are given by 

. 2mHH y 

5 =-—(IeT' ( 34 ) 



mHH x 

7 = 7^T' (35) 

where m is the dimensionless constant that parameterizes the amplitude of the deviation from uj( de ' and can be 

given real values. 

Now, subtracting equation (13"2"j) from (1551) . we get 



Using equations ©, flU]), ([34]) and <(35j) , from (J36J) we obtain 



B 
B 



3(n 2 -l) + m(n + 2) 2 



3(n-l) 
The general solution of equation (|3"T)l has the form 



|J=0, (37) 



B= (fcit + fco)^^ 11 , (38) 

where m = 3iVi + m(n + 2) 2 with Ni = (n- l)(n + 2). 
For A in this case we find 

A = (ht + ko)^ 11 . (39) 

It is important to note here that we obtain power law solution by assuming proportionality relation between 
shear scalar(er) and expansion (9) which seems to describe the dynamics of Universe from big bang to present 
epoch while a series of works: Yadav and Yadav (2011); Yadav et al (2011); Amirhashchi et al (2011); Kumar 
and Yadav (2011); Akarsu and Kilinc (2010) and recently Yadav (2011) have obtained the power law solution 
by assuming special law of variation of Hubble's parameter. So, we represent the new features of power law 
expansion. In this paper, we show how a oc model with metric ([1} behaves in presence of perfect fluid and 
anisotropic DE components. 
Now, the metric (p} reduces to 

ds 2 = -dt 2 + (kit + /c )'""i dx 2 + (kit + k Q )^~ (dy 2 + dz 2 ) . (40) 

In view of the assumption lu 1 -™ 1 ' — const., equation (fi~5j) can be integrated to obtain 

p(">) = p a- 3 ^ (m)+ V t (41) 

where po is the positive constant of integration. 

The physical parameter such as directional Hubble parameter (H x ,H y or H z ), average Hubble parameter (H), 
anisotropy parameter (A m ), shear scalar (a), expansion scalar (6) and spatial volume (V) of model (HO]) are 
respectively given by 

3»(»-l) fa 
m (kit + k ) 

rr - rr - 3 ( n ~ 1 ) kl an 

H v- H '~ ni (M + fco)' (43) 

H=^ (h ** , (44) 

m (kit + ko) K ' 



2(n-l) 2 

(n + 2f 



A m = ;_ , J 2 , (45) 



a«-4fe^, (46) 



n 2 (M + M 2 ' 
3JVi h 



(47) 



m (kit + ko) 

V = (k 1 t + k ) (3Nl/ni) . (48) 
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Figure 1: Plot of matter density (p*-™ 1 ') and DE density (p( de ^) versus time. 



The average scale factor (a) and DP (q) are found to be 

a={k 1 t + k ) iNl/ni) , 

l = 2 + -^(n + 2) 2 



(49) 
(50) 



- ni 

q ~ Nx " ' Ni 

As one sees, the DP q is a constant. The sign of q indicates whether the model inflates or not. A positive sign of 
q, i. e. n\/N\ > 1 corresponds to standard decelerating model whereas negative sign of q, i. e. < n\/Ni < 1 
indicates acceleration. The recent observations SN la, reveal that the present Universe is accelerating and the 
value of DP lies somewhere in the range —1 < q < 0. It follows that in the derived model, one can choose the 
value of DP consistent with observations. 



From (|41|) we than find the energy density of perfect fluid 

p(™)=p (fc li +fc )-3(- < '" ) +l)iV 1 /n 1; 

on account of (|5ip from J7J we obtain the dark energy density as 



Jde 



2 1.2 



) _ 9{2n+l)(n-l) 2 k 



A) 



n 2 (M + fco) 2 (fcii + fco) 3 ^"^ 1 )^/" 1 ' 

Now from (fM|) and (|3"5")l skewness parameter are obtained as 

6m(n-l)iVifc 2 



5{t) = 
7(«) 



9(2n + l)(n - l) 2 fc 2 - p n 2 (fcit + fc )2-3( w c»)+i)jv 1 /„ 1 
3mn(n — l)A^ifc 2 



9(2n + l)(n - 1) 2 & 2 - pon?(M + fc )2-3(«("O+i)Wi/»i ' 
The EoS parameter of DE is given by 

3(n-l) 2 fc 2 (6-2m(n+2)-rii/(n-l)) 



(de 



W 



It is observed that at t 



(feli + fco) 3 ^*™^ 1 '™!/"! 



nfiktt+ko) 2 



9(2rc+l)(ra-l) 2 fc 2 
n^fcit+feo) 2 



(fc 1 t+fc ) 3 (- m + 1 ) Jv i/'-' 



(51) 
(52) 

(53) 
(54) 

(55) 



-Tp, the spatial volume vanishes while all other parameters diverge. Thus the 
derived model starts expanding with big bang singularity a.t t — —j 1 which can be shifted to t = by choosing 
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Figure 2: Plot of EoS parameter of DE components (o/ rfe )) versus time 



ko = 0. This singularity is point type because the directional scale factors A(t) and (B(t)) vanish at initial 
moment. From Fig. 1, we observe that p*-™) as well as p( de ) remains positive during the cosmic evolution. 
Therefore the weak energy condition (WEC) as well as null energy condition (NEC) are obeyed in the derived 
model. Further p^ and p( de ) decrease with time, and approach to a small positive values at the present epoch. 
The parameter H, a and 9 start off with extremely large values and continue to decrease with expansion of 
universe. Fig. 2 clearly shows that a/ de ) evolves with negative and it's range is in nice agreement with large 
scale structure data (Komatsu et al 2009). 
The density parameters of perfect fluid and DE are as follows: 



OiVi /Ci 



n {de) = 



3(2n+l) 



p n 



i2(fcii + fco) 2 ~ 3(a/m)+1)JVl/ni 



(n + 2) 2 3N?kl 
Adding equation (|5B)) and (J3T)) . the overall density parameter (17) is obtained as 



2(n-l) 2 



A, 



(56) 

(57) 

(58) 



This shows that the overall density parameter (17) depends on the anisotropy parameter (A m ). Fig. 3 demon- 
strates the behavior of density parameters in the evolution of Universe with appropriate choice of constants of 
integration and other physical parameters using reasonably well known situations. We observe that initially the 
ordinary matter density dominates the Universe. But later on, the DE density dominates the evolution which 
is probably responsible for the accelerated expansion of present-day Universe. 



4 The statefinder &; Distance Modulus Curves 

Salmi et al (2003) proposed a cosmological diagnostic pair {r, s} called state finder, which is defined as 
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Figure 3: Plot of density parameters versus time. 



The dynamics of statefmder {r, s} depends on constant n\ and N\. It follows that in derived model, one can 
choose the pair of statefmder which can successfully differentiate between a wide variety of DE models including 
cosmological constant, quintessence, phantom, quintom, the chaplygin gas, braneworld models and interacting 
DE models. For example if we put n\ = 0, the statefmder pair will be {1, 0} which yields the AC DM (cosmolog- 
ical constant cold dark matter) model. The statefmder diagnosis for holographic DE model in non flat Universe 
has been analyzed by Setare et al (2007). 



The distance modulus is given by 



\x = 5 log cIl + 25 , 



where the luminosity distance d^ is defined as 



ri(l + z)a 



where z and ao represent red shift parameter and present scale factor respectively. 
Let us now assume that T — kit + k$. Thus equation (141?)) may be rewritten as 

a = T n2 , 



(61) 
(62) 

(63) 



where ri2 = N\/n\ 

For determination of r\ , we assume that a photon emitted by a source with co-ordinate r — r± and T = T\ and 

received at a time Tq by an observer located at r = 0. Then we determine r\ from 



T\ 



T 



Tt 



clT 



(64) 



Solving equations (|6ip — (|64[) . one can easily obtain the expression for distance modulus (p) in term of red shift 
parameter (z) as 

n 2 fci 



fi = 5 log 



H (l - n 2 )(l + z)~ 



(u + ^r 



i 



+ 25. 



(65) 



The comparison between the derived model and SNLS type la supernovae data can be seen in Fig. 4. The dotted 
line represents the observed distance modulus by SNLS type la supernovae data where as solid line represents 
the analyzed distance modulus /i of the derived model. It is observed that the derived model is best fit with 
high redshift values. 
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Figure 4: Distance modulus as a function of the redshift to the derived model compared with SNLS type la 
supernovae data from Astier et al. (2006) 



5 Conclusion 

In this paper, we have studied a spatially homogeneous and anisotropic LRS Bianchi-I space time filled with 
perfect fluid and anisotropic DE possessing dynamical energy density. Studying the interaction between the 
ordinary matter and DE will open up the possibility of detecting DE. It should be pointed out that evidence 
was recently provided by Abell-Clustcr A586 in support of interaction between DE and dark matter (Bcrtolami 
et al 2007; Le Delliou et al 2007). Let us now concentrate on some other works on minimally interacting perfect 
fluid and DE models of Universe, especially the works by Akarsu and Kilinc (2010) and Yadav (2011). In both 
the works, the field equations have been solved by using special law of variation of Hubble's parameter which 
yields the constant value of DP where as the present investigation is one with power law solution by taking into 
account the proportionality relation between shear scalar (a) and expansion scalar (9). It is to be noted that our 
procedure of solving the field equations is altogether different from what Akarsu and Kilinc (2010) have adapted 
in LRS Bianchi-I space-time. However, a common feature of all these power law solution models is to describe 
the dynamics of Universe from big bang to present epoch. 

In the derived model, the EoS parameter of DE (o/ de )) is obtained as time varying and it is evolving with 
negative sign which may be attributed to the current accelerated expansion of Universe. Also note that the 
isotropic distribution of DE is not possible in LRS Bianchi Type I space-time because the isotropic distribution 
of DE leads to the positive value of DP which can not explain the current accelerated expansion of Universe 
while for anisotropic distribution of DE, DP evolves with negative sign. The distance modulus curve of derived 
model is in good agreement with SNLS type la supernovae for high redshift value which in turn imply that the 
derived model is physically realistic. 



The age of Universe is given by 

ni ° Jfei 



Tq = ~ H o 



which is different from the present estimate i.e. To = Hq 1 — HGyr. But if we take fco = and Ni/m — 1, i.e., 
m = 2(1 — n)/(2 + n), where m is the constant describing the anisotropy of DE and n is the constant giving 
the proportionality condition between shear and expansion scalar, then the derived model is in good agreement 
with the present age of Universe. 
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